Current induced entanglement of nuclear spins in quantum dots 
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We propose an entanglement mechanism of nuclear spins in quantum dots driven by the electric 
current. The current accompanied by the spin flip in quantum dots gradually increases components 
of larger total spin of nuclei. This entangled state drastically enhances the spin relaxation rate 
of electrons, which can be detected by measuring a leakage current in the spin-blockade region. 
This mechanism is not relevant in optical experiments which examine the spin relaxation of single 
excitations. 
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The spin relaxation time of electrons in quantum dots 
is an important issue. The time needs to be sufBciently 
long for the implementation of quantum computing de- 
vices utilizing electron spins in quantum dots Q. The 
relaxation time has been found to be quite long by op- 
tical H, ^ and transport experiments In the for- 
mer experiments, the transient measurements indicate a 
quenching of the spin relaxation on the exciton lifetime 
scale in self-assembled quantum dots. In the latter exper- 
iments, the tunneling current is measured through single 
or double quantum dots in "spin blockade" situations, 
where electrons cannot be transported unless the spins 
are flipped in the dot. The relaxation time is as long 
as 200 /IS from a spin-triplet excited state to spin-singlet 
ground state in single quantum dots In weakly cou- 
pled double quantum dots, the current suppression has 
been observed when two electrons occupy the lowest en- 
ergy level in each dot with parallel spins: an electron 
tunneling from one dot to the other is forbidden by the 
Pauli exclusion principle [1). 

In the present paper, we theoretically study a small 
leakage current in spin-blockade regions. The current 
is accompanied by the spin flip in the quantum dots. 
As a spin-flip mechanism, we investigate the hyperfine 
interaction between electrons and nuclear spins [|[ 
Recent experimental results have implied its important 
roles in the leakage current We examine a forma- 
tion of entangled states of nuclear spins driven by the 
electron transport. This mechanism is analogous to the 
current induced dynamic nuclear polarization (DNP) in 
quantum Hall systems |l^ . The DNP is created by the 
electron scattering between spin-polarized regions, which 
leads to a hysterisis in the longitudinal resistance ||l^, ^ . 
In our case, the entanglement of nuclear spins drasti- 
cally enhances the spin relaxation rate of electrons in 
quantum dots, which could be observed experimentally. 
This mechanism is not relevant in the optical experiments 
which examine the spin relaxation of single excitations. 
Hence the role of the hyperfine interaction in the spin re- 
laxation could be quite different in transport and optical 
experiments. 

Model: To consider the spin blockade in quantum dots, 
we adopt a following model. A quantum dot is weakly 



coupled to external leads L,R through tunnel barriers. 
The Coulomb blockade restricts the number of electrons 
in the dot to be A^oi or N^x + 1. Nc\ electrons form a 
background of spin singlet. From lead L on the source 
side, an extra electron tunnels into the dot and occupies 
a single level with envelope wavefunction The spin 

of the electron is either ] {S^ = 1/2) or I {S^ = -1/2) 
with equal probability (we assume an easy-axis of elec- 
tron spins in z direction for a while). The electron stays 
for a long time by the spin blockade |jl^ because the spin 
relaxation time is much longer than the tunneling time. 
After the spin flip in the dot, the electron immediately 
tunnels out to lead R on the drain side, and then the 
next electron is injected from lead L. 

An electron occupying orbital 'il'{r) interacts with N 
nuclear spins, Ifc, by the hyperfine contact interaction. 
N ^ 10^ in GaAs quantum dots. We assume nuclear 
spins of 1/2 for simplicity. The Hamiltonian in the dot 
reads 

N 

i/M = A^^;o|V(rfe)|'S-Ifc«2aS-I, (1) 

fc=i 

where vq is the volume of the crystal cell and I — X^^i ^fc 
is the total spin of nuclei. We have assumed that V'(r) 
is independent of the nuclear site, r^. This is a good 
approximation for a large part of nuclear spins in the 
quantum dot, except in the vicinity of nodes of V'(r), 
since the distance between nuclei is much smaller than 
the size of the dot, or an extension of ip{r). 

Basic idea: The Hamiltonian (|^) indicates that N nu- 
clear spins interact with a common "field" of electron 
spin although there is no direct interaction between them 
(dipole-dipole interaction between nuclear spins is weak 
and neglected). This field results in the entanglement 
among the nuclear spins. To illustrate this, let us con- 
sider the simplest case of = 2 and begin with a po- 
larized state of nuclear spins \Iz — \/2)i\Iz — 1/2)2- 
An electron with spin | [) tunnels into the dot and is 
spin-flipped by the hyperfine interaction. Then the state 
of nuclear spins becomes (| — l/2)i|l/2)2 + |l/2)i| — 
l/2)2)/V2 = |J = 1,M = 0), where J and M are 
the total spin and its z component, respectively. This 
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is an entangled state. After the electron tunnels out 
of the dot, the next electron is injected with | |) (or 
I t)). The spin- flip probability of the second electron 
is proportional to |(l,Tl;T (Dl-^hf |1, 0; i (T))P = ^a" . 
This value is twice the probability in the case of non- 
entangled states |l/2)i| - 1/2)2 or | - l/2)i|l/2)2. In 
general, the capability of state \J,M) to flip an electron 
spin is (X (J ± M){J + 1). 

In presence of easy-axis: We formulate the entangle- 
ment of N nuclear spins. In the initial state, they are 
not entangled at all and oriented randomly 

- (ci|l/2>i+di|-l/2>i)®... 

® (cAr|l/2)jv + rfjv|-l/2)jv). (2) 

{cfc,c?fc} are randomly distributed (|cfcp -I- — 1). By 
the transformation of the basis set from {|/z)i ® \Iz)2 ® 



• • • } to the total spin of all the nuclei, {| J, M)} ( J < N/2, 
~J<M<J), Eq. (1) is rewritten as 

^(0)^ E (^J%i,x\J^M,\), (3) 

.J,M,X 

where the coefficients {Cj^l^ ^} are randomly distributed 

JM A \^JM A I ~ Index A distinguishes states with 
the same J and AI . The number of such states is given 

by 

iV' 

KiN, J) = (2 J + l) ^ ^+y+2 ),(^-2j -jT- (4) 

After the injection of the first electron, say, | I), the time 
evolution of the dot state by is 



-tajt/h 



,J,M,\ 



2J+ 1 



J - M + 1 + (J + M)e'"(2J+i)t/Rl I A) ® I i) 



iajt/h 



+ E ^fkxViJ + M){J -M+ l)^y--(l - e-(2^+i)*/'')| J, M - 1, A) | T) 



J,M,X 



2J+1 

r 



(5) 



The expansion to the lowest order in at/h yields the spin- 
flip probability, P^") = F^^^at/hf, with 

^(0)^ E \cfM,x\\J±M){JTM + l) (6) 

J,M,A 

(the lower sign indicates the electron spin flip from 
I t> to I i>). When > 1, \Cfl.iJ^ can be re- 
placed by 1/2^ (law of large number). Then F^^^ — 
{1/2'')J:j,mK{N,J){J ± M){J T M + 1) = A^/2. 

= {N/2){at/hf is identical to the spin-flip prob- 
ability which would be evaluated on the assumption that 
one of the nuclear spins is flipped with the electron spin. 

The electron tunnels off the dot immediately after the 
spin is flipped. Then the wavefunction ^''"^ shrinks to 

x|J,Mt1,A). (7) 

is partly entangled since it contains more compo- 
nents of larger J and smaller \M\. The degree of the 
entanglement increases each time an electron is injected, 
spin-flipped, and ejected out of the dot. After n events. 



the state of nuclear spins becomes 

^in) ^ ^ c(J jj,M,A), (8) 

J,M,A 

pin) ^ ^ \d;%^^\^j±M){jTM + i){io) 

J,M,X 

F*^"^ is expressed as 

/'"^ = ^EW^)[(J±M)(JTM + l)r(ll) 

J,M 

Figure 1 shows the distribution of (a) the total spin J 
and (b) its z component M. The spin-flip processes in- 
crease the weight of larger J and smaller \M\. Although 
M changes less faster than J, this means that the total 
nuclear spins tend to develop in the plane perpendicular 
to the easy-axis of electron spins. 

This entangled state enhances the spin-flip probability. 
For the [n + l)th electron, the probability is given by 
pin) ^ F^'^^at/hf. P(") with flxed t is shown by sohd 
line in Fig. 2(a). Using Eq. (|l|), we flnd that F^") « 
(A^/2)n for K n < A^/2, and F^") w {N/2f for A^/2 < 
n. Therefore the probability increases with n linearly 
(p(n) nP(°'^) and flnally saturates (P(") « (A^/2)p(o)). 
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FIG. 1: The distribution of (a) the total spin J and (b) its 
z component, in the state of nuclear spins. An easy-axis of 
electron spins is assumed, (a) p(J) = X^m a I'-'jAf aI^i ^'^'^^ (b) 

p{M) = X] J A I'-^jA/ aP' where n is the number of transported 
electrons accompanied by the spin flip. The number of nuclear 
spins is A'^ = 100. In Eq. we take the geometrical mean 
between upper and lower signs. 
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FIG. 2: (a) The spin-flip probability P'") with fixed t, as a 
function of n (number of transported electrons accompanied 
by the spin flip), (b) The spin-flip current I{t) as a function 
of time t. The cases in the presence and absence of easy- 
axis of electron spins are indicated by solid and broken lines, 
respectively. The number of nuclear spins is N = 100. In the 
presence of easy-axis, we take the geometrical mean between 
upper and lower signs in Eq. (]ll|). 



In absence of easy-axis: Until now, we have as- 
sumed the presence of easy-axis of electron spins. In 
the absence of the axis, the spin of the nth elec- 
tron is oriented in an arbitrary direction, cos 0'"-' | t 

) + sin0(")e*'^ I i). In this case, the previous formu- 
lation can be applied after the rotation of {C'j"2/a}' 
jC/(J,0("),(^("))M,A./'Cg,„A' to align 



the z axis parallel to the spin direction. If the spin direc- 
tion of incident electrons is completely random, we can 
average over the ^-component of the total spin. Then 
the coefficients {Cj"^^ in Eq. are replaced by the 

averaged values {C'f'l}. They develop by 



C 



(") 

J, A 
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J(J+ 1) 
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(n-1) 
J, A 



(12) 



where 



}^Y.K{N,J){2J+l) 



'-J{J+l] 



.(13) 



As in the previous case, components of larger J increase 
with n, which enlarge the spin-flip probability by the 
factor of 2J(J+l)/3. The probability is given by P(n) = 
G''"^(at//i)^, which is shown by broken line in Fig. 2(a). 
This behaves as P(") « (2n/3)p(") for 1 < n < N/2, 
and « (Af/3)F(°) for A^/2 < n, with fixed t. 

Leakage current: The enhancement of the spin-flip 
probability is reflected by the leakage current. Instead 
of calculating the current in the previous situations, now 
we consider the case in the presence of magnetic field: 



Zeeman energy for electrons is much larger than the hy- 
perfine interaction, whereas that for nuclear spins is neg- 
ligible. Then iJhf can be treated as a perturbation. The 
magnetic field makes an easy-axis of electron spins. We 
also take into account the electron-phonon interaction for 
the energy conservation at the spin flip The Hamil- 
tonian in the dot reads H = Hd + iJph + -ffhf + ^^oi-ph, 
where Hc\ = -i?z(a|aT ~a|aj, i/ph = J2q ^q^q^qi and 

being creation and annihilation operators for an electron 
in the dot (phonon) . The spin-flip processes in the lowest 
order in iJhf + ^^oi-ph (flrst orders in and in Hd-ph) 
change the nuclear spin state in the same way as in Eq. 

® 0- 

The spin-flip rate for the (n -I- l)th electron is T^"-' cx 
^("), xhe time interval between nth and {n + l)th elec- 
tron transports is A<^"^ « l/r^'^\ Thus the current is 
given by I{t) = eT^") at t = Y.]Lo ^^^^^ • This current 
is shown in Fig. 2(b), as a function of t (solid line). The 
approximate form of F^"^ yields an asymptotic form of 
the current 



lit) 

m 



e(A^/2)rW 



[t < tsat) 
(isat < t). 



(14) 



The current grows with time drastically and finally 
saturates. The saturation time is given by tsat — 
ln(iV/2)/r(0), where F^") is the spin-flip rate of an elec- 
tron with the non-entangled state of nuclear spins. The 
increase in the current continues until the entangled state 
is broken by dephasing effects. 

Even in the absence of magnetic field, the electron en- 
ergies for I t) and | |) could be split in Hd, e.g. by 
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spin-singlet or -triplet formation with an electron trapped 
in the other dot, in double dot systems. In this case, 
the easy-axis of electrons would be absent (if the crys- 
tal field is weak enough). The perturbation with respect 
to ffhf + -ffci-ph changes the nuclear spin state following 
Eq. (|T|). The spin-flip rate is r(") cx G^"^ The leakage 
current I{t) is represented by broken line in Fig. 2(b). It 
is written approximately as 



m 
m 



e(2r(o)/3)e2r 
e(A^/3)r(0) 



(0) 



t/3 



[t < isat) 
(isat < t), 



(15) 



where i^at = 31n(7V/2)/(2r(°)). The enhancement of the 
current is less prominent than in the presence of easy- 
axis. 

Discussion: We have proposed the current induced en- 
tanglement of nuclear spins in quantum dots. The cur- 
rent accompanied by the spin flip in the dots gradually 
increases components of larger total spin of nuclei, which 
significantly enhances the spin relaxation of electrons. As 
a result, the leakage current in the spin-blockade region 
grows drastically with time and finally saturates. This 
mechanism is not relevant in the optical experiments 
which examine the spin relaxation of single excitations. 

The saturation time of the leakage current is of the 
same order as l/F^'^-', spin relaxation time with non- 
entangled state of nuclear spins, which is ~ 100 ns or 
much larger |Q, The current is enhanced during the 
dephasing time T2. When the entangled state is broken, 
the current is suppressed. Possibly this would result in 
a current fluctuation with time which has been observed 
recently The characteristic time of the fluctuation is 

A possible origin of the dephasing is the dipole-dipole 
interaction between nuclear spins. Note that, among 
N nuclear spins participating in the formation of total 
spin, I J, M), this interaction conserves the total spin, and 
hence does not destroy the entangled state. The interac- 
tion between one of the N spins and a spin surrounding 
the dot results in the dephasing. Estimating T2 in our 
case and analyzing evolution of nuclear spins after the 
dephasing are beyond the present paper. 

We discuss the validity of our simple models. We have 



disregarded the spatial variation of the envelope wave- 
function of electrons in the quantum dot. Although the 
entanglement of nuclear spins is generally seen in the 
presence of the variation (see Eq. (4) in Ref. [^), the 
enhancement of the spin relaxation is more effective with 
larger total spin J. Hence the capability to flip an elec- 
tron spin is determined by the effective number Nos of 
nuclear spins which feel the same electron field, ip{r) = 
const., since Jmax = Ncs/2. The generalized evolution of 
nuclear spins is an interesting problem. 

If the contribution from higher energy levels of elec- 
trons in quantum dots is not negligible, through electron- 
phonon interaction Q or spin-orbit (SO) interaction 
p6[ , it should be taken into account carefully. Partic- 
ularly, the coexistence of hyperfine and SO interactions 
would complicate the evolution of the nuclear spin state 
because the terms of | |) and | t) are mixed on the right 
side of Eq. (||). We have also disregarded higher-order 



tunneling processes 1 17 , which could play a role in the 
spin relaxation of electrons However, the processes do 
not influence the entanglement of nuclear spins discussed 
here. 

Finally, we comment an analogy between our mecha- 
nism and the Dicke effect of superradiance [|l^ |l^ . The 
spontaneous emission of photons is enhanced from N 
atoms with two levels (pseudo-spin Sz = ±1/2) if all 
of them are excited initially. This is due to the forma- 
tion of pseudo-spin state \J,M) with J = A''/2. (Starting 
from I J, J) , the state of N atoms changes like a cascade, 
I J, J), I J, J — 1), • • • , emitting photons.) A similar effect 
has been proposed for the emission of phonons from N 
equivalent quantum dots [|l9|, |2^. The atoms (quantum 
dots) correspond to the nuclear spins in our model, while 
the emission of photons (phonons) to the spin flip of elec- 
trons. A main difference is the initialization. N excited 
states have to be prepared by the pumping in the Dicke 
effect, whereas the initialization is not necessary in our 
mechanism. 
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